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ABSTRACT 

We present a ^-band spectroscopic study of the stellar and gas kinematics, mass distribution, and stellar 
populations of the archetypical starburst M82. We used the ^^CO stellar absorption bandhead at 2.29 fim 
(CO2.29) to measure the rotation curve out to nearly 4 kpc radius on both the eastern and western sides of the 
galaxy. In contrast with the nearly Keplerian gas dynamics suggested by previous measurements of H I and CO 
emission from the ISM, which imply a truncation in M82's dark matter halo, our data show that the rotation 
curve is in fact flat on 1 - 4 kpc scales. The kinematics of the Br7, H2, and He I emission lines are consistent 
with, although characterized by slightly higher velocities than, the stellar kinematics. We derived M82's mass 
distribution from our stellar kinematic measurements and estimate its total dynamical mass is ~10^^ Mq. We 
measured the equivalent width of CO2.29 (W2.29) as a function of distance from the center of the galaxy to 
investigate the spatial extent of the red supergiant (RSG) population. The variation in W2.29 with radius clearly 
shows that RSGs dominate the light inside 500 pc radius. M82's famous superwind is likely launched from this 
region, where we estimate the enclosed mass is < 2x 10^ Mq. 

Subject headings: galaxies: individual (M82) - galaxies: starburst - galaxies: kinematics and dynamics - 
infrared: galaxies 



1. INTRODUCTION 

Observations suggest that star formation in galaxies is regu- 
lated by feedback processes. In rapidly star-forming galaxies, 
feedback operates on large scales, driving galaxy- wide super- 
winds, which can eject the raw materials for future stellar gen- 
erations if the interstellar material is launched with enough 
momentum to escape the host galaxy's gravitational poten- 
tial. Otherwise, the ejected material may be re-accreted in a 
so-called "fountain flow". These issues are crucial for under- 
standing the formation of stars in rapidly star-forming galax- 
ies across cosmic time. 

The exceptionally high star formation rates present in star- 
burst galaxies make them optimal targets to probe feedback 
processes. The focus of this paper is the archetype of this 
class of object s, M82. This galaxy ' s proximity (3.63 Mpc ; 
V = 17.6pc; Freedma netal.( 11994 iSakai & Mado?a [T999I: 



Gerke et al. 2011) and nearly edge-on geometry (/ = 80°; 
McKeith et al., 1993) make it a superb laboratory for studying 



the physics of galactic winds and star formatio n. At infrared 
wavelengths, its luminosity is 6 x 10^^ Lq (ISanders et al.l 
|2Q03|) , which makes it one of the brightest and most observed 
infrared objects in the sky. Most of its infrared luminosity 
originates in the starburst core, which is se verely obscured by 
dust at visible and ultraviole t wavelengths (iRieke et al.ll 19801: 
iForster Schreiber et al.ll2QQli) . 

The origin of M82's starburst is thought to be a tidal in- 
teraction with its massive companion, M81, ^ 10^ years ago 
iGottesman & WeHachew"1977'; ' Lo et al.|[T987l: iTelesco et all 
[l991; Achtermann & Lacy 1995). The pronounced HI brid ge 
that connects the galaxies (lCottrelllll977l: lYun et al.|[T993h is 
strong evidence for such an interaction. In addition to the 
dynamical effects associated with tidal forces from M8 1 , the 
global dynamics and evolution of the starburst are further per- 
turbed by the presence of a stellar bar (Telesco et al. 1991), 
which may have formed as a result of the interaction. 

As an interacting, barred galaxy, M82's observed dynami- 
cal properties are distinct from those seen in typical Sb and Sc 



galaxies. Molecular gas measurements by 'Young & Scovill3 
(1984) show that the CO velocity field contains warps and pe- 
culiarities that are consistent with a significant non-circular 
velocity component, perhaps due to radial infall or a polar 
ring tilted with respect to the major axis. In addition, mea- 
surements of CO (2,1) bv lSofue et al.1 (11992 ) suggest that the 
rotation of the gas disk is nearly Keplerian at radii > 1 kpc. 
Sofue (1998) reproduced this Keplerian motion with a model 
in which the outer disk and halo were stripped during the en- 
counter with M8 1 . Although flat rotation curves are observed 
in the vast majority of galaxies, Rubin & Ford (.1982) found 
falling rotation curves in several tidally interacting galaxies, 
which suggests that some dark matter is stripped from galax- 
ies during interactions. 
IForster Schreiber et al] 



(I2003h proposed that M82's star- 
burst activity occurred in two successive episodes that each 
lasted a few million years and peaked about 10^ and 5x 10^ 
years ago. In this model, the encounter with M8 1 triggered the 
first burst in the central 500 pc by inducing strong large-scale 
torques, loss of angular momentum, and an infall of molec- 
ular clouds towards the nuclear region. The second burst 
occurred in a circumnuclear ring and along the stellar bar, 
and is attributed to bar-induc ed dynamical resonances. The 
IForster Schreiber et al.l (|2003|) starburst model can explain the 
apparent abundance of red supergiants (RSGs) that dominate 
the light within 500 pc radius. 

In order to understand the dynamics of the ejected gas 
and constrain the physics of feedback processes in M82, 
it is essential to accurately measure the total dynamical 
mass on both the small scales of star-forming regions and 
the large scales of the entire galaxy. Rotation curves are 
the majo r tool for d etermining the distribution of mass in 
galaxies (|Sofue & Rubin 2001), and many studies have em- 
ployed gas and stellar kinematics at a ran ge of wavelengths 
to measure radial velocities in M82 (e.g Saito et ajj [T9841: 
Goetzetal. 1990; McKeith etal. 1993; Achtermann & Lacvl 
IT995uSofue et al.i 19921) . iGreva (i2011i) includes a compilation 
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of these results. Given that M82 has recently undergone tidal 
interactions and harbors an extraordinary superwind, stellar 
kinematics should be a more reliable tracer of the mass dis- 
tribution than the gas. Nearly all previous stellar kinematic 
measurements used visible- wavelength tracers, which may be 
compromised by the substantial dust attenuation present in the 
galaxy. The near-infrared /sT-band is a better tracer of the stel- 
lar kinematics, as attenuation in the ^-band is approximately 
a factor of ten less than in the V-band. In this study, we use 
^-band spectroscopy to study M82's gas and stellar kinemat- 
ics, mass distribution, and stellar populations. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained ^-band spectra of M82 at the Mount Graham 
International Observatory with the LUCI-1 spectrograph and 
the Large Binocular Telescope (LBT) on 2011 February 10. 
The observations were made through a ^-band filter (2.05 - 
2.37 /im) with the N1.8 camera, which has a plate scale of 
0'''.25 pixel"^ and a field of view of 4' x 4\ A grating with 210 
lines mm~^ was used with a slit width of V\ This combination 
yields a spectral resolution of A/ A A ^ 3,000 or 100 km s~^ 
Two slit positions, which we designate as M82E and M82W, 
were used to measure the rotation curve out to 4'(~4kpc) 
from the center of the galaxy. We centered each end of the slit 
on the nuclear region at a position angle of 67.5° (Figure [TJ. 
The observations were made with 300 second exposures in an 
object-sky sequence for a total on-source integration time of 
1 hour per slit position. A telluric standard close in air mass 
to M82 was observed before and after both the M82E and 
M82W observations. 

We carried out standard data reduction within the 
IRAlQ software package. The two-dimensional data were 
transformed onto an orthogonal wavelength/slit-position grid 
based on a wavelength solution from the sky's bright OH 
emission lines. This transformation linearized the dispersion 
to ^1 .64 A pixel"^. We used the xtellcor (IVacca et al.l 
l2003h package within IDL to correct for absorption from the 
Earth's atmosphere. The telluric standard spectra had slightly 
higher spectral resolution than was obtained for M82. This 
is likely because the light from the stars did not fill the 1" 
wide slit. This complication was remedied by convolving the 
telluric spectra with Gaussians to match their resolutions to 
M82. 

3. ROTATION CURVE 

To measure the rotation curve, we symmetrically extracted 
thirty apertures each from the M82E and M82W slit positions. 
We progressively increased the bin sizes at larger radii to 
achieve a sufficient signal-to-noise ratio in the stellar contin- 
uum. Figure |2] shows illustrative ^-band spectra from M82's 
nucleus and eastern side. We measured the rotation curve out 
to nearly 4 kpc radius with the ^^CO stellar absorption band- 
head at 2.29 /im (CO2.29) and to ~1 kpc radius with gas emis- 
sion (He I, H2, Br7) from the interstellar medium. Because of 
the extreme tidal forces and galactic wind in M82, the gas and 
stars may trace distinct kinematic structures, and we discuss 
each in turn. 

3.1. Stellar Absorption 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 
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Fig. 1 . — Digital Sky Survey image of M82 with overlays of our two slit 
positions. Each slit is 4^ long and they overlap on the center of the galaxy. 
The widths of the slits have been increased by a factor of six for ease of 
visibility in this figure. 



Q 



M82E 

He I H2 


19" (-330 pc) radius 




M82E 


8.8" (-150 pc) radius 

Bry 


Nucleus 


12CO 



2.0 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 



A [/jm] 

Fig. 2. — ^-band spectra from M82's nucleus, 8.8'' east from center, 
and 19'' east from center. CO2.29, which is the dominant absorption feature, 
originates in the atmospheres of cool, evolved stars. The He I and Br7 lines 
trace star formation, and the H2 is evidence of shocked gas. The strength of 
the emission lines peaks near 19". 

CO2.29 originates in the atmospheres of cool, evolved stars 
and provides a strong absorption feature for stellar veloc- 
ity measurements (e.g. 'Gaffn ey et al.l 119951) . Previous stud- 
ies of M82 have used Ha and Na D to measure the rotation 
curve with stellar kinematics (e.g. lMayalll[T96Ql: iGoetz et aP 
|199Q|) . yet these measurements are likely compromised by 
the Ay ^ 5 ^ 25 attenuation toward the central region 
(iRieke et al.lll98Ql :lLester et al.'199Q: Puxlev'T991). At some- 
what longer wavelengths, McKeith et al. (1993) measured the 
stellar kinematics with Ca II stellar absorption lines (8662, 
8542, 8498 A). If the near-infrared extinction wavelength de- 
pendence is given by the relation Ax (x A~^^^ (lDraindll989h . 
attenuation for CO2.29 is approximately a factor of five less 
than Ca II. With the copious amount of gas and dust present, 
an improvement at this level may reveal kinematic features 
that are otherwise obscured. 

To measure the stellar rotation curve, we selected six aper- 
tures with high signal-to-noise as cross-correlation templates 
for relative velocity measurements. We used the rvs ao pack- 
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Fig. 3. — Stellar rotation curve based on C02.29- The horizontal error bars indicate the aperture size of the measurement, and the vertical error bars represent 
the uncertainty of the velocity measurement. All data points without error bars have uncertainties less than 2 km s~^ and aperture sizes less than The flatness 
of the rotation curve is inconsistent with Keplerian motion. The symmetric, sharp features on either side of the nucleus are due to the non-circular motions of the 
stellar bar At ~1 kpc on either side of the galaxy, there is an apparent "velocity reversal", which may be due to symmetric spiral arms. 



age within IRAF to cross-correlate all sixty apertures against 
each of the six templates. We restricted the wavelength range 
from 2.28 to 2.365/im to focus the analysis on CO2.29. This re- 
sulted in six independent rotation curves, all of which we nor- 
malized to be symmetric around zero velocity. After correct- 
ing for M82's nearly edge-on inclination, we took the mean 
velocity at each point to be the value of the rotation curve 
at that point, and the uncertainty was calculated as the rms 
between the six curves. Figure [3] shows our stellar rotation 
curve. The horizontal bars represent aperture size. All data 
points without error bars have uncertainties less than 2 km s~^ 
and aperture sizes less than 5''. Table[T]lists the measurements 
(apertures selected as cross-correlation templates are marked 
with a dagger). 

The observed stellar velocities show the stellar rotation 
curve is flat on scales of 1 - 4 kpc. This is in stark contrast 
with the nearly Keplerian gas dynamics sugges ted by H I an d 
CO emission from the interstellar medium (So fue et al.ll 1992b . 
The gas, however, is severely disrupted by the wind and may 
also be affected by radial inflow and tides (lYoung & Scovilld 
[1984). The stellar data are therefore expected to trace the 
mass distribution much more reliably. 

The rotation curve for M82E appears to decline slightly 
from ~ 2 ^ 4 kpc. This decline has ~2cr significance and 
may be due to tidal stripping caused by the interaction with 
M8 1 . The declining rotation curve is at the approximate lo- 
cation of a flaring in the stellar disk and the northern H I 
streamer (Yun et al.lll993l) . No decline is seen in the M82W 
data. Measurements at larger radii are needed to determine 
w hether or not the fall off is due to tidal stripping. 

iTelesco et al.l (1199 Ih discovered M82's stellar bar with ob- 
servations at infrared wavelengths and estimated it to be ap- 



proximately 1 kpc in length. Numerical simulations have 
shown that galaxy interactions are effective at inducing galac- 
tic bars (e.g. Noguchi 1988; Miwa & Noguchi 1998), and the 
strong evidence for past tidal interaction observed in M82 
makes it an excellent testing ground for theories of bar for- 
mation. Although it is difficult to observationally determine 
the origin of M82's bar, the effects of its presence can be seen 
in our stellar rotation curve as the symmetric, sharp features 
on either side of the nucleus. 

The second bump at ~500 pc may be the result of a super 
star cluster (seen in the two-dimensional data) that fell on the 
slit . This cluster, which was identified as J0955548+694953 
by llMcCrady et al.l (l2QQ3h . likely dominates the flux from this 
aperture, and its velocity excess with respect to the main rota- 
tion curve likely indicates strong non-circular motion. How- 
ever, it is puzzling that multiple points define the feature and 
that no other studies have seen this second peak in the rota- 
tion curve. Other possible explanations for the feature are the 
presence of tightly wound spiral arms and patchy ^-band ex- 
tinction. 

Our data contain evidence of a "velocity reversal" (veloc- 
ity decrease) at ^ 1.0 kpc on b oth the eastern and western 
sides of the galaxy. iGrevd (1201 1|) highlights the existence of a 
100 km s~^ velocity reversal of gas and/or stars at ~ 1.0 kpc in 
US and He stellar absorption on the eastern side of the galaxy, 
as well as in ionized gas emission lines at visible and blue 
wavelengths on the eastern and western sides; interestingly, 
the feature is apparently absent in near-infrared Ca II stel- 
lar absorption. Although the velocity amplitude of our de- 
tection is approximately a factor of ten less, the reversal is 
present in our data. Local extinction is not likely to be the 
culprit, as the feature is symmetric with respect to the nu- 
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TABLE 1 
Stellar Rotation Curve Data 
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All values are measured with respect to M82's center. The positions are flux weighted. At M82's distance, l' 



i 1.056 kpc. 



T Cross-correlation aperture 

cleus. iGrevd (|2Q11|) investigates the hypothesis that the re- 
versal is due to l ocal surface density enhancements associated 
with spiral arms (Mayya et al. 2005), but finds inconsistencies 
with this picture. For instance, the absence of the velocity re- 
versal in near-infrared Ca II stellar absorption suggests that 
it is not inherent to the disk, and the blue- wavelength detec- 
tions of the eastern velocity reversal reach negative velocities, 
which cannot be explained from the rotation of a stellar disk 
with embedded corotating spiral arms. Nevertheless, the de- 
tection of the reversal in ^^CO stellar absorption and the sym- 
metry about the nucleus indicate the feature is present in the 
stellar population and implies it is due to a local density en- 
hancement, such as spiral arms. 

3.2. Gas Emission 

There are several prominent gas emission lines in our data 
that are characteristic of starburst activity, most notably Br7, 
He I, and H2 (Figure [J]). The Br7 and He I are recombina- 
tion lines that trace photoionized nebulae, and the H2 is a 
ro vibrational line that traces shocked gas. In contrast to the 
observed strength of CO2.29 relative to the continuum, which 
decreases with distance from the nucleus, the emission lines 
show an increase of intensity with distance inside 350 pc ra- 
dius and a rapid decrease with distance beyond 1 kpc radius. 
This is consistent with the i ntensity variations measured by 
iForster Schreiber et al.1 (I200L) . 

For each aperture listed in Table [T] with sufficient signal- 
to-noise, we used the emsao package within IRAF to mea- 
sure the velocity offset of each emission line relative to its 



rest wavelength (ICoxII200Q|) . We normalized the emission line 
rotation curves to the CO2.29 velocity in an aperture ^1.25^' 
west of the nucleus. Figure |4] compares our gas and stellar 
rotation curves with the nearly Keplerian gas dynamics mea- 
sured in CO (2,1) by lSofud (Il998i) . The emission lines used 
to measure the gas kinematics are only detected out to 1 kpc. 

On average, our measurements show that the gas moves 
with higher velocities than the stars inside 1 kpc. This can 
be explained by forces from the stellar bar, w hich t end to 
drive the gas towards the center of the galaxy. Sofud (Il998|) 
measured even larger velocities with radio observations of 
CO (2,1) and HI in this region Measurement s of opt ical O II 
emission from ionized gas by iMcKeith et al.l (Il993h . on the 
other hand, yielded lower velocities than the stars. As the 
depth of an observation is a function of wavelength, these dis- 
crepancies may be the result of optical depth variations. 

4. MASS DISTRIBUTION 

The mass of M82 has a much larger contribution from 
atomic and molecular ga s than seen in typical spiral galaxies. 
lYoung & Scovilld (119841) estimate that between 30% and 40% 
of the galaxy's mass is in the ISM and that the H2/HI mass ra- 
tio is approximately 10:1. These authors postulate that the 
high gas fraction may be due to infalling material that con- 
tributes to the available reservoir. The estimated dynamical 
mass of M82 ranges from several times 10^ 10^^ MQ_jn 
the literature (e.g . Burbidge et al. 1964; Crutcheret al. 1 9781: 
'Young &Scoviriail984: iGoetz etal.iil99Q: iSofue et al.lll992h ; 
however, the substantial dust attenuation and non-circular mo- 
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Fig. 4. — The gas and stellar rotation curves from this work with the nearly Keplerian measurements of CO (2,1) bv lSofud (119981) . The bottom panel is a 
magnification of the dashed region in the top panel. On average, the gas has larger velocities than the stars. The stellar rotation curve is not consistent with 
Keplerian motion. 

tions associated with tidal forces and the stellar bar make the 
estimation of the mass challenging and uncertain. 

Figure [5] shows the mass profile based on our stellar kine- 
matic measurements. Similar to the studies mentioned above, 
we assumed circular motion to calculate the enclosed mass 
at each point. While non-circular motions are expected in 
the central region, and in particular in the bar region, circu- 
lar motion is a good assumption at large scales. At 4 kpc the 
mass distribution appears to be asymptotically approaching 
~10^^ Mq. This is likely the total mass of the galaxy, as these 
measurements are based on stellar motions and extend out to 
nearly 4 kpc from the nucleus. 

Circular motion is not a good assumption in the cen- 
tral region, where the bar and perhaps other asymmetries 
in the potential produce non-circular motions. Neverthe- 
less, here we estimate the mass within 500 pc is ^ 2 x 
10^ M0, which is consistent with mass measurements o f this 
region by Youns & Scovilld (11984 . iGoetz et al.l (119901) . and 
iGreve et al. ( 2002). This mass estimate may be high by a 
factor of two based on the mass distribution for M82E and 
M82W shown in Figure [5] 

5. STELLAR POPULATIONS 

The CO absorption index in stars increases with increas- 
ing stellar luminosity, decreasing effectiv e temperature, and 
increasing metallicity ("Frogel et al.lll978L il98 3). Large CO 
indices are indicative of either a stellar population domi- 
nated by RSGs or one dominated by old, metal-rich gi- 
ants. These two possibilities have led to much debate on 
the nature of M82's evolved stellar popula tion, particularly 
in the nucleus (iForster Schreiber et al.ll200l]) . Starburst mod- 
els with RSG-dominated nuclei, howeve r, provide the best 
fit to the observed properties of M82 (iRieke et al.l 119931: 
IForster Schreiber et al. .20031 




1.5 2.0 2.5 

Radius [kpc] 

Fig. 5. — Mass profile derived from our stellar rotation curve. Circular 
motion was assumed in the derivation. The horizontal error bars denote the 
aperture size of the measurement. The apparent decrease in mass at short radii 
may be due to the non-circular motions associated with the stellar bar. The 
superwind is likely launched from the region inside the dashed line, where 
we estimate the mass is < 2 x 10^ Mq. 

Our data reveal very deep CO absorption in M82's starburst 
core. In Figure [6l we show the equivalent width of CO2.29 
(W2.29) as a function of radius. The data clearly show a sig- 
nificant increase in W2.29 within 5 00 pc radius. Based on 
the diagnostic dia grams proposed bylOriglia et al.l (P^93) and 
lOliva et al] (Il995|) . these results suggest that RSGs dominate 
the near-infrared c ontinuum throughout th e starb urst core. At 
yet smaller radii, IForster Schreiber et al.l (12001 ) found that 
W2 29 within the central 35 pc is characteristic of RSGs with 
r,// = 3600 -4200 K. 
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The dramatic contrast in W2.29 between the starburst core 
and the stellar disk can be expl ained by the succ e ssive star- 
burst model proposed by Forst er Schreiber et alJ (l2QQ3h , in 
which the disk and the starburst core are composed of two 
distinct stellar populations. 
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Fig. 6. — W2.29 as a function of distance from the center of the galaxy. 
The horizontal error bars indicate the aperture size of the measurement. The 
peak in W2.29 inside the dashed region is indicative of a stellar population 
dominated by RSGs. 



6. SUMMARY 

We have used longslit, ^-band spectroscopy from the LBT 
to measure the stellar and gas kinematics and study the mass 
and stellar population distribution in M82. The primary points 



of this paper can be summarized as follows. 

We used ^^CO stellar absorption to measure the rotation 
curve out to nearly 4 kpc on the eastern and western sides of 
the galaxy. While previous measurements of HI and CO from 
the interstellar medium suggested that M82's rotation curve is 
Keplerian, our data show that it is flat on scales of 1 - 4 kpc 
(Figure O. The rotation curve also contains the signatures of 
the stellar bar and an apparent "velocity reversal" at ~ 1.0 kpc 
on either side of the galaxy that may be the manifestation of 
symmetric spiral arms. 

We measured the rotation curve with the Br7, He I, and 
H2 emission lines from the interstellar medium out to 1 kpc 
from the nucleus. These gas rotation curves have velocities 
10-50 km s~^ greater than the stellar velocities (Figure |4l). 
The gas dynamics are expected to be disrupted more than the 
stars by the presence of tides and galactic winds, and may thus 
be a less accurate tracer of the total mass distribution. 

We used the stellar absorption data to measure the mass dis- 
tribution (Figure [5]) and investigate the stellar populations as a 
function of radius (Figure[6]). We estimate the total dynamical 
mass of M82 is ~ 10^^ Mq. The strong variation in W2.29 with 
radius clearly indicates that the starburst is located within the 
central 500 pc radius. The superwind is likely launched from 
this region, and we estimate its mass is < 2 x 10^ Mq. This 
mass estimate is somewhat uncertain due to the presence of 
non-circular motions, but nevertheless provides an upper limit 
to the gravitational potential of the starburst region. 

We are grateful to Jeff Blackburne and Jill Gerke for taking 
the data as part of the OSU/RC queue. JPG is grateful for 
support from an undergraduate research scholarship from the 
College of Arts and Sciences at the Ohio State University. PM 
is grateful for support from the NSF via award AST-0705170. 
TAT is supported in part by NASA grant NNXIOADOIG and 
an Alfred P. Sloan fellowship. 
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